Quantum Destruction of Spiral Order in Two Dimensional Frustrated Magnets 
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We study the fate of spin— 1/2 spiral ordered two dimensional quantum antiferromagnets which 
are disordered by quantum fluctuations. A crucial role is played by the topological point defects 
of the spiral phase which are known to have a Z2 character. Previous works established that a 
non-trivial quantum spin liquid phase results when the spiral is disordered without proliferating the 
Z2 vortices. Here we show that when the spiral is disordered by proliferating and condensing these 
vortices, valence bond ordering occurs due to quantum Berry phase effects. We develop a general 
theory for this latter phase transition and apply it to a lattice model. This transition potentially 
provides a new example of a Landau-forbidden deconfined quantum critical point. 
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I. INTRODUCTION 

In recent years we have learnt a great deal about 
possible quantum phase transitions out of Neel ordered 
phases of spin- 1/2 quantum antiferromagnets in two spa- 
tial dimensions^. Much of the work have focussed on 
circumstances where the Neel ordering is collinear (i.e., 
where all the spins ahgn along a common axis, but, al- 
ternate in direction). It is now known that such collinear 
Neel states can give way to paramagnetic valence bond 
solids (VBS), which break lattice translation symmetry 
but not spin rotation symmetry, through generic sec- 
ond order quantum phase transitions. Since the two 
phases (Neel and VBS) break incompatible symmetries of 
the Hamiltonian, a continuous transition between them, 
within the conventional paradigm of critical phenom- 
ena (known as the Landau-Ginzburg- Wilson paradigm' ), 
is not allowed except at fine-tuned multi-critical points. 
Hence the above generic second order transition has been 
dubbed as Landau-forbidden^. The notion of deconfined 
criticality has been introduced in the context of this 
and other similar phase transitions"". Possible continu- 
ous phase transitions between the collinear Neel phases 
and gaplcss spin liquids that do not break any symme- 
tries have also been described"^. 

In contrast, quantum phase transitions out of the spi- 
ral ordered Neel phase (i.e., with non-coUincar spin pat- 
terns), that may occur on non-bipartite lattices, have not 
been as extensively studied. Early theoretical works'*'"' 
discussed continuous phase transitions from such a spi- 
ral state to a gapped Z2 spin liquid with fractionalizcd 
spin-1/2 spinon (bosonic) excitations and associated non- 
magnetic Z2 vortices (dubbed visons*^). These studies 
have given rise to speculations that disordering spiral 
magnets naturally lead to deconfined spin liquids and 
the transition between the spiral phases and conventional 
paramagnets with valence bond order has remained un- 
explored. This is despite good theoretical motivations 
and some numerical evidence that exists for understand- 
ing such transitions as is reviewed below. 



Quite generally one can view the nature of the tran- 
sition out of the spiral from the perspective of its topo- 
logical defects. These are the well-known point-like Z2 
vortices in two spatial dimensions^. If the spiral is de- 
stroyed without proliferating these vortices a gapped 
Z2 spin liquid with fractionalizcd bosonic spinon excita- 
tions emerges and the resulting transition is rather well 
understood- the transition is described in terms of the 
condensation of the bosonic spinous, while the visons re- 
main gapped throughout; the critical point belongs to 
0(4) universality class''. 

However, the more conventional transition out of the 
spiral is driven by proliferating the Z2 vortices, such that 
both the spinous and the visons become critical at the 
transition. What is the nature of the resulting param- 
agnet and the associated transition? This is the funda- 
mental question that we will formulate an answer to in 
this paper. We will argue that the quantum Berry phase 
effects lead to VBS order in such paramagnets and the 
associated transition may be generically second ordered. 
A schematic picture for a quantum phase transition be- 
tween the spiral and the dimer is shown in Fig. 1. In this 
paper we shall outline the general structure of the field 
theory for a generic continuous transition between the 
two phases and work out an example in context of a lat- 
tice model. We emphasize that the theory developed here 
can capture similar transitions on other lattices. Exam- 
ples of such extensions are considered towards the end 
of this paper. As pointed out earlier, since the spiral 
and the VBS break different incompatible symmetries of 
the Hamiltonian, a generic second order transition be- 
tween them is an example of violation of the conventional 
Landau-Ginzburg- Wilson' theory phase transitions. 

The organization of the rest of the paper is as follows. 
We begin, in Sec. II, by briefly reviewing some of the 
motivations that exists for studying quantum phase tran- 
sitions out of spiral magnets. There are a number of nu- 
merical studies of frustrated spin-1/2 quantum magnets 
on a square lattice. Most pertinent to this paper are 
the exact diagonalization calculations on the "Ji — J3" 
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model on a square lattice®. Over a range of J3/J1 (see 
below) the system orders into a spiral pattern. In a dif- 
ferent, though proximate, range a paramagnetic ground 
state apparently obtains. This later state is believed to 
have VBS order. Similarly, as we argue below, the same 
model on a rectangular lattice possibly has a direct quan- 
tum phase transition between the spiral and the VBS. 

To construct the critical theory, we find it convenient 
to start from the spiral phase and consider disordering it 
by proliferating its defects. To this end we discuss differ- 
ent aspects of the spiral order parameter in Sec. Ill-the 
order parameter manifold and its topological defects, the 
Z2 vortices. It is known that such vortices carry non- 
trivial Berry phases. The role of the vortices in the tran- 
sition is effectively captured using the well-known spinon 
parametrization of the spiral order parameter. This is in- 
troduced in Sec. Ill A. The lattice field theory consistent 
with different symmetries of the problem is then written 
down in Sec. IV and different limits of this theory are dis- 
cussed. It is shown that in presence of the Berry phase, 
the paramagnet obtained by proliferation of the Zi vor- 
tices is indeed dimerized. The critical theory, containing 
both soft spinon and vison modes, is discussed in Sec. V. 
The spinous and the visons see each other as source of 
TT-flux (mutual semions) . This arbitrarily long range sta- 
tistical interaction between them is implemented through 
a mutual Ising Chern-Simmons interaction term. After a 
series of transformations, the final theory is written down 
in terms the spinous coupled to a U(l) gauge field. The 
gauge field is compact and the VBS order parameter is 
proportional to the monopole operator of the gauge field. 
The critical theory contains doubled-monopole operators 
where the gauge flux changes by ±47r. However, combin- 
ing available analytical and numerical results, we argue 
that this doubled-monopoles are dangerously irrelevant 
right at the critical point, i.e., while they are relevant in 
the paramagnetic phase, they are irrelevant at the criti- 
cal point. Thus the gauge field is rendered non-compact 
at the critical point. This field theory, expected to de- 
scribe the critical point, is an anisotropic version of the 
NCCP3 field theories (see below). The issue of exten- 
sion of the present theory to other lattices is briefly taken 
up in Sec. VA. Finally we conclude with our remarks in 
Sec. VI. The details of various calculations are summa- 
rized in Appendices A-D. 



II. LATTICE MODELS 

We are interested in spin- 1/2 frustrated quantum anti- 
ferromagnets in two spatial dimensions described generic 
Hamiltonians 
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where S.^ is the spin- 1/2 operator at the ith lattice site 
and J > denotes nearest neighbour antiferromagnetic 



FIG. 1. Schematics of a generic continuous quantum phase 
transition between the spiral and the valence bond paramag- 
net. The transition is brought about by the proliferation and 
condensation of the Z2 vortices (refer to the text) achieved 
by tuning the parameter A. 



exchanges. The ellipses denote other spin rotation invari- 
ant interactions that may be tuned to drive the phase 
transition between the spin ordered and the paramag- 
netic phases. Though our formulation is gencrically ap- 
plicable to all cases of direct and continuous quantum 
phase transition between the spiral and the dimer, here, 
for concreteness, we consider a simple lattice model to 
illustrate our results. Extensions of these ideas to other 
relevant lattice systems are quite straight forward and 
indicated towards the end of the paper. 

Consider spins on a rectangular lattice with first ( Ji) 
and third neighbour (J3) antiferromagnetic exchanges 

Uf = ( Ji Sr • Sr+x + Sr • Sr_|_2x) 
r 

+ A ( Ji Sr • Sr+y + J3 Sr • Sr+2y) (2) 

r 

where r denotes the sites of the rectangular lattice. The 
couplings along y direction arc A (anisotropy factor) 
times those along the x direction (0 < A < 1). The 
two tuneable parameters in the Hamiltonian are A and 
{•hi Ji)- The presence of the latter leads to frustration 
of exchange interactions. 

Various limits of this model are well known. For A = 
0, one has decoupled spin chains with nearest and next 
nearest neighbour couplings, where, it is known that for 
Js/ Ji < 0.241 there is power-law Neel order, while above 
this true long range VBS order is obtained^. On the 
other hand, A = 1 represents an isotropic square lattice 
with nearest and third nearest neighbour interactions. 
Numerical results'^ '^'^ suggest that this has at least three 
phases. For J3/J1 < 0.3 the usual collinear Neel state 
obtains while for J3/ Ji > 0.7 the ground state shows non- 
coUinear spiral order. At intermediate values one gets a 
paramagnet which possibly breaks lattice symmetry, i.e. 
a VBS state. 

Imagine sitting in the spiral phase at A = 1 (by choos- 
ing J3/ Ji > 0.7). What happens if A is decreased towards 
zero? Clearly decreasing A increases quantum fluctua- 
tions so that the spiral order will be destroyed below some 
critical Ac. For very small A the VBS order of the decou- 
pled chains will persist as columnar dimer order with a 
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2-fold degenerate ground state. Could this VBS state 
persist all the way up to Ac ? If so could the resulting 
transition be second order? While the first question can 
only be answered by detailed numerical studies of this 
particular microscopic model in future we will formulate 
an answer to the second one in this paper. 

Having established the setting we now look at the spi- 
ral phase and consider disordering it through prolifera- 
tion and condensation of its defects. 



III. THE SPIRAL ORDER PARAMETER AND 
ITS TOPOLOGICAL DEFECT 

Consider a magnetically ordered phase where the av- 
erage spin orientations are coplanar. Such a spin pattern 
is characterized by 

(Sr) - (ni(r) cos (Q • r) + n2(r) sin (Q • r)) ^ 0, (3) 

where, Q is the ordering vector which may be commensu- 
rate or incommensurate with the underlying lattice struc- 
ture. For Q 7^ nn in being an integer) the spin order is 
non-collinear and the forms a spiral pattern (as shown in 
Fig. 1). ni(r), n2(r) are two mutually orthogonal unit 
vectors, i.e., 

na{r) ■ np{r) ^ Sap. (4) 

Here a, (3 = 1,2. The main difference between the 
collinear and spiral order parameter is that in the former 
only a single unit vector (ni or is sufficient, while the 
latter requires both rii and n2. The spiral order param- 
eter is then given by an §0(3) matrix 

3? = [ni, n2, ni x . (5) 

This difference in the order parameter manifold has 
several important consequences, the foremost being three 
spin-wave modes in the spiral phase compared to two in 
the collinear phase. Another feature, which is central 
to this paper, is the fact that topological defects of the 
spiral are distinct from that of the collinear Neel order. 
For the spiral phase, the order parameter has the sym- 
metry of three-dimensional rotation group SO(3). The 
order parameter manifold is isomorphic to /Z2, which 
is an unit sphere in 4-dimensional space whose antipodal 
points are identified. Thus the order parameter space is 
compact and doubly connected}^ . 

In two spatial dimensions, this allows topologically sta- 
ble point vortices that are characterized by a Z2 quantum 
number^ 

Hi (53/^2) =^2, 

U2 {S^/Z2) = 0, (6) 

where Hi and 112 denotes the first and second homotopy 
groups respectively^''. These Z2 vortices are distinct from 
the 0(2) point vortices in a familiar two-dimensional 



Xy-model. In case of the Z2 vortices, they can have 
a winding number which is either or 1 (i.e. the winding 
number is defined modulo 2). This is a refiection of the 
fact that there are only two classes of closed path in the 

/ Z2 manifold and all closed paths can grouped in one 
of these two classes. Hence, the Z2 vortices are charac- 
terized by an Ising quantum number. Two such vortices 
can annihilate each other. 

In the ordered phase the energy of a single vortex di- 
verges logarithmically with the system size and free vor- 
tices are absent. However the vortices play a crucial role 
when the spiral order is destroyed by condensing them. 
To this end it will be extremely useful to set up an ef- 
fective description of the spiral and proximate param- 
agnetic phases that captures easily the role of the Z2 
vortices. The formulation in terms of ni_2 along with the 
constraints (Eq. 4) is not convenient for this purpose. 
Instead it is useful to introduce an alternative parame- 
terization of the spiral order parameter. 



A. The spinon parametrization of the spiral order 
pEirameter and the Ising gauge fields 

We introduce the well-known redundant description of 
the order parameter in terms of the spinon variables. 
Specifically we write^'^"* 

n+ = ni + m-2 = ea^zpa-a-yZ^, (7) 

where crs are the three Pauli matrices and is the 
two dimensional completely antisymmetric matrix (ei2 = 
-£21 = -l,eii = £22 = 0). z = izi,Z2) is the two 
component complex spinon field with unit norm, i.e., 

Z^Za = 1- (8) 

The spinous arc thus bosons that transform as spin-1/2 
under spin rotations. We note that the spinous defined 
here are related to the Caley-Klein parameters (which in 
turn are related to the Euler angles) used in describing 
the motion of a classical rigid body^^. 

At this point it is useful to consider the global symme- 
tries of the effective action written in terms of the spinon 
fields. The action must be invariant under global SU(2) 
spin rotation under which the spinous transform as 

Za VapZp, (9) 

where Vap is an §U(2) matrix. In addition there are sev- 
eral other discrete symmetries of the Hamiltonian- vari- 
ous lattice transformations and time reversal. For a rect- 
angular lattice, the transformations of the spinon under 
these symmetries are given by : 

T ■ r ^ p-jQ-a/2 

T:za^ iz*a, (10) 
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Where Ta represents unit translation along a, Jig denotes 
reflection about an axis of symmetry of the lattice, 3 de- 
notes inversion of the lattice and T denotes time reversal. 
The Translation symmetry imphes that the continuum 
theory has an extra U(l) symmetry along with the spin 
SU(2) and other discrete symmetries. 

The spinon parameterization is two-to-one and there 
is a discrete Z2 gauge symmetry corresponding to the 
change of sign of the z fields independently at each site 



z(r) — z(r) 



(11) 



which leaves the order parameter invariant. This reit- 
erates the fact, in terms of the spinous, that the order 
parameter manifold is S'^/Z2. 

Around a vortex the order parameter fields Hq, are sin- 
gle valued. However the spinon fields, Za, changes sign. 
Thus the spinon wave function changes its sign on going 
around a vortex. In other words the spinous and the Z2 
vortices see each other as sources of tt flux. 

Following Lammert et al}"^ ^ a fruitful description of the 
vortices is now achieved by introducing an Ising gauge 



field cr,: 



±1, minimally coupled to the spinous, on 



the links of the direct lattice. The Z2-vortices are then 
associated with the magnetic flux of a Z2 gauge field 



(12) 



where the product is taken over the links of the plaquette 
of the Ising gauge field. = — 1(-|-1) indicates the pres- 
ence(absence) of a vison inside the plaquette. It is im- 
portant to note^^ that visons are well defined excitations 
even in the paramagnetic phase. Here, the spinon fields 
fluctuate wildly and so does their corresponding paths in 
the order parameter space. However, for spinous describ- 
ing a closed loop around a vison, the paths necessarily 
end at diametrically opposite points in order parameter 
manifold (S'V^2). 

The Berry Phase of the Z2 vortices 

The above classical picture must be augmented with 
the correct quantum Berry phase term. Semi-classical 
analysis^^'^^ for the spiral phase shows that the non- 
trivial Berry phases are solely associated with the vor- 
tices and are given by 



exp 



E (1 

i,j—i-\-r 



n < 



(13) 



where (Ji^i+r are the Z2-gauge fields on the time-like links 
of the (2 -I- 1)D space-time lattice. A different study^°, 
starting from the spin disordered phase, recovers the 
same Berry phase term. Eq. 13 is also, not coinciden- 
tally, the Z2-Polyakov loop term obtained in the analysis 
of the quantum dimer models^^ . 



It is useful to introduce^^ a new set of Ising gauge field 
/i^^* = ±1 (not to be confused with the Hab introduced 
later. We have used same notations as in Ref. 16 for easy 
reference) on the links of [ah) dual (2-1-1) dimensional 
lattice such that 



n 



ext 
t^ab 



-1 



(14) 



for all the space-like plaquettes of the dual lattice. Then 
one can write the Berry phase in Eq. 13 as 

ITT \r-^, 
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(i-^..)(i-n/^a?)' 



(15) 



□ 



where the links (ij) and (ab) belong to the direct and the 
dual lattices respectively such that the dual plaquette is 
pierced by the link (ij). 



IV. LATTICE FIELD THEORY 

Various global symmetries and the Z2 gauge structure 
now fully determine the effective action which is invariant 
under the transformation group 

[SU(2) X U(l)]^,„,,, X [Z2]„ (16) 

plus discrete lattice symmetries and time reversal (see 
Eq. 10). The minimal imaginary time Ginzburg-Landau 
action in (2 -I- 1) space-time lattice, consistent with these 
symmetries, is 

S = S.+§s, (17) 

where 

§z = (7ij (^zj ■ Zj + h.cj — r (^zl ■ Zj — zt ■ z. 



(18) 

and §B is the Berry phase contribution given by Eq . 13. 
Eq. 17 may be derived more formally using Hubbard- 
Stronovich transformation as shown in Appendix A. In- 
tegrating out the higher energy spinous generate several 
terms allowed by the symmetry, the foremost being the 
Z2-Maxwell term which gives kinetic energy to the gauge 
fields 



p 



(19) 



where the sum is taken over the plaquettes of the space- 
time lattice and Kp depends on the plaquette orienta- 
tion. Thus the effective low energy action is given by 

= + §B, (20) 

where Sz, So- and $b are given by Eqs. 13-19. This action 
given by Eq. 20 is powerful enough to capture the large 
part of the phase diagram including the spiral, dimer and 
the Z2 spin liquid phases and the associated transitions. 
However, in this paper we shall focus exclusively of the 
possibility of the direct second order transition between 
the spiral and the dimer phases. 
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The VBS Phase 



The spiral phase is obtained by condensing the spinons, 
i.e. (zq) 7^ (hence (n) ^ 0). On the other hand, 
deep inside the paramagnet the spinon excitations are 
gapped and can be integrated out from Eq. 20 to yield 
the effective action 



S' 



(21) 



This is the action for the Odd Ising gauge theory (Odd- 
IGT)-^^ which have been studied extensively in context of 
quantum dimer models^^ . In (2+1) dimensions, Odd-IGT 
is fruitfully studied using Ising duality^''''^'' connecting 
the Odd-IGT to the fully frustrated transverse field Ising 
model on the dual lattice^^. The dual action is given by 



3< 



FFTFIM 



J^KatPlt^lf Pi ^TY, Pi- 



rn) 



Here are dual Ising spins that occupy the sites of the 
dual lattice. The Berry phase (Eq. 13) imposes the con- 
dition of full frustration in this dual model through /i^^* 
as defined in Eq. 15. We note again that these fields 
obey the constraint Yla Mab* ~ ~^ '^^'^^ space-like dual 
plaquettes^*^. Without loss of generality we take the dual 
exchange couplings Kab > 0. 

As derived in Appendix B, the vison creation operator, 
V'Lonl"). is proportional to pj, i.e., 



V'Lon(a) ~ Pa- 



(23) 



Wc note that the vison resides in the plaquettes of the 
direct lattice which is equivalent to the sites (a) of the 
dual lattice. Also the VBS order parameter, '^vBsi'ij) is 
proportional to {PaPb)^ 



VBS 



(pIpD. 



(24) 



where, as described in Appendix B, the link (ij) on the 
direct lattice crosses the link (ab) on the dual lattice. We 
note that The VBS order parameter is bilinear in vison 
operators and is gauge invariant as one expects. 

The dual Ising spins undergo an ordering transition 
and the ordered phase is characterized by 



(pl) + 0. 



(25) 



However due to the frustration, the ordered phase breaks 
lattice translation symmetry^*'. Eq. 25 indicates that the 
visons proliferate and condense in this phase. The VBS 
order parameter gains a non-zero expectation value 



4- 



VBS 



(26) 



and reflects the broken lattice translation symmetry. 
Thus the vison condensate indeed describes a dimerized 
paramagnet^^. 



FIG. 2. The Fully frustrated model on the rectangular lat- 
tice: The black bonds (both solid and dashed) are ferromag- 
netic and the blue (wavy) bonds are antiferromagnetic. The 
gauge is chosen in such a way that around each plaquette that 
rin Mab* = — 1 is satisfied. Further the magnitude of the solid 
black bonds are taken to be 77 while those of the dashed black 
bonds and the blue bonds are taken to be 1, where 77 > 1. 
The unit-cell, as shown in the figure, contains two sites. 



The critical vison modes 

The momenta of the vison modes that condenses de- 
scribes the VBS ordering pattern near the transition. 
Following Blanckstein et al.^^, these modes can be ob- 
tained from a soft mode analysis of Eq. 22. In our case, 
the dual Ising model is described on the dual rectangular 
lattice. There is a gauge redundancy in choosing p.^* 
in Eq. 22 and maintaining the constraint Yia PaT ~ ~^ 
over all space-like plaquettes (see discussion following Eq. 
22). In the rest of this calculation we shall assume, with- 
out any loss of generality, a specific gauge to identify the 
soft modes'^^. This is shown in Fig. 2. There are two 
soft modes: Aq and A^ at wave-vectors (0,0) and (7r,0) 
respectively. 



Ao 



1 



A^ 



(27) 



(Note that the modes arc un-normalized. Also see Fig. 
2 for definition of the dual anisotropy factor rj.) 

The low energy vison modes near the transition are 
linear combination of Aq and A^ : 



$(r, t) = «'o(r, <)Ao + *^(r, t)A„, 



(28) 



where 5'o(r,i) and ^'7r(r,i) are the two complex ampli- 
tudes. The effective Ginzburg-Landau action for the soft 
vison modes can be constructed by considering transfor- 
mations of the two amplitudes under various symmetries 
of the Hamiltonian. Due to the gauge redundancy in 
choosing /i^^*, the transformations we need to consider 
are the projective symmetry transformations (PSG) of 
the soft vison modes^^. It is useful to define a single 
complex vison mode 



= 







(29) 
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The Projective symmetry transformations for the soft vi- 
son modes are given by: 



Jy ^ -I**, 



(30) 



Here Tx and Ty are unit translation afong x and y respec- 
tively, 3?7r is rotation by tt about a lattice point, Rx,Ry 
are reflections about x and y axes respectively and 3 is 
Lattice Inversion. The lowest order action allowed by 
these transformations is given by 



{a.b) 



fc)-r„2jcos(40a). (31) 



The modes transform as "XY spin" under different sym- 
metries of the Hamiltonian up to 4*'' order. As outlined 
before, in terms of the original spins, the ordered phase 
of this theory represents the VBS and e"^ is identified 
with the vison creation operator. 

Having outlined the critical spinon and the vison 
modes we now proceed to construct the critical field the- 
ory. 



V. THE CRITICAL FIELD THEORY 

At the critical point both visons and spinous are gap- 
less. Further, they see each other as sources of tt flux 
(mutual semions). Hence, there is an arbitrarily long 
range statistical interaction between them. Such interac- 
tions may be effectively implemented by introducing two 
Ising gauge fields crij(/iah), on the links of the direct (dual) 
lattice, coupled by an Ising Chern-Simmons term^^ : 



Scs 



ITT 
T 



(ab) □ 



(32) 



Here we note that this mutual Ising Chern-Simmons cou- 
pling is different from the mutual U{\) Chern-Simmons 
term used for similar problems in Ref. 30. We also take 
this opportunity to distinguish between the approach of 
the present paper and that of Ref. 30, which is based 
on the mutual U{1) Chern-Simmons theory. In particu- 
lar, we focus on the latter's limitation in capturing the 
possible direct second order transition between the spiral 
and the dimer. In their work Xu et al. introduced the 
critical spinous and the visons like this work. However 
in their continuum theory they use two U{1) gauge fields 
coupled via a Chern-Simmons term to implement the mu- 
tual semionic statistics between the spinous and the vi- 
sons. This is valid, as the authors acknowledge, when 
the number of spinon and vison excitations are small'^''. 
However near the spiral to dimer transitions there is a 
proliferation of low energy spinon and soft vison excita- 
tions and the U(\) theory breaks down. On the other 
hand the mutual U{\) Chern-Simmons theory success- 
fully captures the transition between the spiral and a Z2 



spin liquid or collinear Neel phase. Indeed, the phase di- 
agram obtained in that formalism shows that the spiral 
and the dimer phases are always separated, either by a Z2 
spin liquid or a collinear Neel phase. A direct transition 
between the two requires fine tuning to a multi-critical 
point and hence is not generic. On the other hand the 
present formalism can tackle the situation as we outline 
below. 

Following the short digression we now return to the 
discussion of the critical theory. The critical Landau- 
Ginzburg action is 



(33) 



where the different terms are given by Eqs. 18, 31 and 
32 respectively. The long range interaction between the 
spinous and the visons, as encoded in Sc5, makes the 
analysis of this field theory difficult and a series of trans- 
formations are required to cast the theory in a useful 
form. To start with, we neglect the effect of the terms 
with coefficient r (in Eq. 18) and F,, (in Eq. 31) and 
consider their effects later. 

Within Villain approximation^^''^^, the first term in the 
vison action (Eq. 31) can be written as (details of the 
calculations are given in Appendix C) 



If ' 

(ab) 



(34) 



where Jab is the integer valued vison current defined on 
the links of the dual space-time lattice. The constraint 
on Jab, as derived in Appendix C, are 



V- J = 



and 



(35) 



(36) 



The first constraint is the continuity equation for the 
vison current reflects the fact that for F^, — 0, in Eq. 31, 
the vison number is conserved. The second constraint, on 
the other hand is the restatement of the mutual semionic 
statistics between the spinous and the visons. The zero 
divergence constraint is satisfied by defining an integer 
valued vector field a on the links of the direct lattice 
such that 



J = V X a. 

The second constraint then becomes 

We now write 



(37) 
(38) 
(39) 



where bij is an integer field and Sij = 0(1) when Uij 
even(odd). Eq. 38 then gives cTy = 1 — 2sy . 
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The condition on a^j to be an integer may be imple- 
mented by applying a soft potential: 

Vsoft = -g ^ cos {2TTbij) = -5 ^ cTi] COS (Tray ) (40) 



where g > 0. At this stage it is useful to break the 
2 complex spinon fields, zi,Z2 into a 4-component real 
vector field: 



Z2 = 1^3 + ii^i- 



(41) 



Further, rescaling the gauge potential aij — >■ aij/n and 
choosing the transverse gauge 

V • a = (42) 

by defining a scalar field C on the direct lattice such that 

a,y aij + AQj (43) 

we have 

§C = - X! '^'i t*^'^* cos (fly- + (0 - Cj})] 



So far we have neglected the efi^ect of r coupling present 
in the microscopic model (Eq. 18). Without this term 
the symmetry is enlarged from the original microscopic 
§U(2) X U(l) to §0(4). Hence we expect that on taking 
r 7^ in the microscopic model , terms allowed by the mi- 
croscopic symmetry would be generated on coarse grain- 
ing and reduce this enlarged symmetry back to the mi- 
croscopic one. Such a generic term can be found starting 
from the universal covering group §U(4) of uni-modular 
4-component complex vectors. 

The §1LJ(4) is generated by fifteen 4x4 matrices 
{Ta,^J.a,Ta^lb} where 



(50) 



where are the usual Pauli matrices and T"-fJ' implies 
direct product. Out of these, the 6 purely imaginary 
antisymmetric ones that generate §0(4) are 

{'Ty, My; '^xMy; '^yf^X7 '^yf^z: ^zl-^y}- (^-^) 

Since §0(4) = SU(2) x §U(2) these 6 generators can be 
broken into 2 mutually commuting sets: 







■ 




■ 




Ta = 















where v = {vi, 1^2,1^3, 1^4) ■ Summing over (Tj^s and collect- 
ing the relevant coupling terms give the critical action to 
the lowest order 



(52) 



E 



(V,j X fly) 



(45) 



The §U(2) group generated by the first set of generators 
represent spin rotations (see Appendix D). Now we can 
choose a generic term that breaks the second §U(2) to 
U(l). This residual U(l) can then be identified with the 
global U(l) symmetry under lattice translation (refer to 
our previous discussion). This is easily done by consid- 
ering the most generic term at the lowest order allowed 
by various point group symmetries. The transformations 
are given in Appendix D with the result being 



where we have introduced complex 4-component vectors 

Xm = J^»ae"^'. (46) 

The constraint over the x fields being that they are uni- 
modular and parallel, i.e., 

X*aXa = 1 
XaX/3 - XpXa = 0. (47) 

This is implemented by incorporating a soft potential: 

VLft^Ml-ixtfix^f) (48) 

where 770 > 0. The continuum limit for this critical theory 
may now be written using a soft-spin description. 



ix* 



xf 



(53) 



(the term X* ■ Ty • % is forbidden by time reversal). So 
the final continuum field theory is 



d^xdr 



id. 



ia^)x\^+p\xf 



u{\x\'Y 



-m{x*?{x? + i{x*ryxY 



-^(Vxa)^ 



(54) 



d X dr 



{dt.-ia^,)x? +P\X? + u{\x\ 



'VoixlHx? 



(V X a) 



(49) 



Thus the final field theory is described in terms of 4 com- 
ponent complex matter fields coupled to a non-compact 
U(l) gauge field. 

Fate of Doubled Instantons: Now consider the role of 
Ft, (in Eq. 31), which introduces the 4-fold anisotropy. 
For Fu = 0, from Eq. 31, the (vison)number operators 
conjugate to the fields are conserved (Eq. 35). This 
is equivalent to the flux conservation of the U(l) gauge 
field (see Eq. 37). Finite F„ destroys this conserva- 
tion. Remembering that e"^ is a vison creation operator, 
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we see that the Ty term allows the simultaneous appear- 
ance/disappearance of 4 visons or change of gauge flux 
by ±47r, i.e., a doubled instanton operator''^. 

In Eq. 54, the condensation of the spinous lead to 
spiral ordering. Once they condense, the gauge field dy- 
namics is gapped through the Anderson-Higg's mecha- 
nism and the instantons are suppressed. On the other 
hand the instantons are relevant in the paramagnetic 
phase and their condensation lead to VBS order. Thus 
this field theory captures the right limits and suggests 
that there may be a direct transition between the spiral 
(Higgs phase) and the dimer (confined phase) 

This direct transition may be continuous only if the 
doubled instantons are irrelevant at the critical point. 
Presently, accurate estimates of the scaling dimension 
of this doubled instanton operator is missing. While 
these may be obtained numerically in future, here we 
make a crude estimate. Large- A^, RPA treatment of the 
gauge fluctuations'^'', suggest that the scaling dimensions 
of instanton of charge q is proportional to q^N'^^. Re- 
cent numerical studies'^'' on the isotropic NCCP^ model 
find that a single instanton has a scaling dimension of 
0.63. Combining these, we find that the scaling dimen- 
sion of the doublcd((7 = 2) instanton operator (A) is 
A = I X 2^ X (0.63) « 5.04 > 3. Hence this naive esti- 
mate suggests that doubled instantons are irrelevant at 
the critical point and so the U(I) gauge flux is conserved 
right at the critical point. This emerging U(l) symmetry, 
absent in the microscopic model, is typical to deconfined 
quantum critical points. An extensive characterization of 
such critical points is given in Ref. 2. 



A. Extension to Other Lattices 

Having fleshed out the details of the calculation for the 
case of rectangular lattice, we now briefly point out the 
extension of this theory to other two dimensional lattices, 
where similar transitions may arise. To this end, we note 
that apart from stabilizing the respective phases, the in- 
formation about the lattice structure enters the analysis 
in the following way: (1) The point group symmetry of 
the lattice dictates the transformation of the spinon and 
only terms invariant under these symmetry transforma- 
tions are allowed in the action; (2) The vison PSG de- 
pends on the point group symmetry of the dual lattice, 
the latter being related to that of the original lattice. 

For the square lattice, repeating the above calculation 
one finds that there arc 2 soft vison modes and the vison 
action is given by a XY model like Eq. 31, but now with 
an eight- fold symmetry breaking term (cos [80]). Prelim- 
inary analysis suggests that the structure of the theory 
remains unchanged except for the fact that we now have 
quadrupled instanton operator which is connected to the 
fact that on the square lattice the VBS state has 4-fold 
degeneracy. 

Another interesting case is that of the anisotropic tri- 
angular lattice, where a similar transition may occur. 



Starykh et al. showed that if one takes spin- 1/2 chains 
with both nearest neighbour and next nearest neigh- 
bour antiferromagnetic interactions and then couple sev- 
eral such chains antiferromagnetically so as to create an 
anisotropic triangular lattice, then such a microscopic 
Hamiltonian allows Spiral and dimer phases in proxim- 
ity to each other'^''. It is interesting to ask whether this 
system allows a direct transition, and if so, what is the 
critical field theory? We can extend our present anal- 
ysis to answer the question of structure of the critical 
theory. Once again preliminary analysis"^^ suggests that 
there are four vison modes and the minimal vison action 
has 0{2) X 0(2) symmetry. 

Extensions to other lattices may be made similarly. 
Finally, it is also interesting to note that for analogous 
cases in 1 D spin chains the visons will always proliferate 
(since in (1 + 1) D the Z2 gauge theory is always in a 
confining phase) leading to dimerization. This helps us 
to sec why the spin-1/2 Ji — J2 chain is dimerized in the 
regime where classical analysis predicts spiral order. 



VI. CONCLUSION 

In this paper we have outlined the field theory for a 
direct second order quantum phase transition between a 
spiral state and a VBS in context of a concrete spin-1/2 
lattice model in two spatial dimensions. This is poten- 
tially a new example of deconfined quantum criticality. 

The continuum theory given Eq. 54 is an anisotropic 
version of the non-compact CP^ model (NCCP^) and 
belongs to the general family of NCCP^~^ critical the- 
ories [N denotes the number of matter components). 
Anisotropic NCCP^ describes the transition between 
coUinear Neel and VBS phases in easy-plane spin-1/2 2D- 
antiferromagnets^, while anisotropic NCCP^ describes 
the continuous transition between spin-nematic and VBS 
in spin-1 2£'-antiferromagnets'^^. Similar to the these 
other field theories describing the deconfined quantum 
critical points, the present one is also a strongly coupled 
theory in (2+1) dimensions. Detailed analysis of such 
theories require numerical investigation. 
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Appendix A: Derivation of the Effective Action 

Here we give an alternate way of arriving at the ef- 
fective action given by Eq. 18. The lattice action writ- 
ten in terms of the spiral order parameters has a global 



9 



SU(2) X U(l) symmetry. The most general action consis- 
tent with this symmetry is introduced by the following 
partition function: 

Zn = y"pjdniidn2i exp(-Sn-SB) (Al) 

i 

§n = (nii • nij + n2i • n2j) (A2) 

m 

where denotes the contribution of the Berry phase 
term. Using the spinon parametrization as given by Eq. 
7, this becomes 

Zz= JJdZiarf^a'^dZjaP - l)exp(-S^ -Sfi), 

(A3) 



Now following the standard slave particle saddle point 
treatment^*^ we approximate the auxiliary field by an uni- 
form saddle point value = Lg (Vij). The saddle point 
respects the global symmetries SU(2) x U(l), but it docs 
not obey the local Z2 gauge symmetry of the spinous. 
The fluctuations of Lij about the saddle point involves 
both amplitude and the phase fluctuations. It is expected 
that the fluctuations that restore the Z2 gauge symmetry 
of the spinous are energetically cheapest. Thus consid- 
ering this class of fluctuations around the saddle point 
which re-implements this gauge symmetry, we write 

Lij = L^aij, (A7) 

where aij = ±1 is a Z2 gauge field defined on the space- 
time lattice. Using this in Eq. A6 we get back Eq. 18. 



Appendix B: Odd Ising gauge theory, visons and 
VBS order parameter 

We now consider the massive amplitude fluctuations of 
Lij in action given by Eq. A6. About the saddle point 
let us replace the fluctuations (Eq. A7) by: 

Lij — {Lq + lij)(Tij, (Bl) 



= -Jl E \A • - ^ (z^ • Zj - Z] . Zi) 

(«j> (ij) 

(A4) 



Introducing auxiliary field Lij on the links and perform- 
ing the Hubbard-Stratonovich decoupling of the first 4- 
spinon term of the spinon action above, we get 



(A5) 
(A6) 



I 

where lij are the amplitude fluctuations of auxiliary flcld. 
Putting this in Eq. A6 we have, 

- j2E(^r^j-^j -^0' ^^^^ 

Finally integrating out the spinons (we can set J2 = for 
this leading order calculation) we generate several terms 
with the leading order term being : 

LQli,i+7i{0'i,i+,i(7i+7j,i+ri+T<7i+'ri+T,i+T(^i+T,i ), (B3) 

i,V 

where the product of the gauge flelds is over the time- like 
plaquette associated with space-like h^i+ri {v being the 
displacement to the nearest space-like neighbour). This 
coupling between lij and as result to 

{h,i+r]) ^ {<^i,i+r]<^i+7),i+ri+T'^i+r]+T,i+TO'i+T,i) (B4) 

or in the time continuum limit of the Ising gauge theory 
we have: 

~ (B5) 

At this point we note that the spinon parametrization 
of Eq. 7 implies that the VBS order parameter 

(Si.Si+,)^(|zrzi+,|2) (B6) 



J Y\_'^^iadz*^ dLijdL*^ 5{\z^a\'^ - l)exp(-S^ - §b), 

i 

^ E l^'^ l' - E [^^1^1 ■ + ^■'') E {^l ■ - 4 • ^i)' ' 

(ij) {ij) (ij) 
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Then from Eq. A6 we find that z[ • Zi_|_,, ^ ^ij+jj- Thus 
wc have: 

(Si • Si+,,) ^ Ll-^+n (B7) 
Thus to the leading order: 

(Si • Si+.;j) ^ Loh^i+ri- (B8) 
Using Eq. B5 we finally get 

(Si . Si+,) ^ (B9) 

We now use two well-known results from the duality be- 
tween Ising gauge theory and Ising model in (2 -I- 1) di- 
mensions (we have taken the time continuum limit as 
before). These are: 

= PaPt (BIO) 

and 

^□=114==/'- (Bll) 

□ 

where, in Eq. BIO, the link (ij) is on the direct lattice 
and (a, b) arc sites on the dual lattice such that the link 
(ab) cross (ij). arc the dual Ising spins as defined 
below Eq. 22 in the main paper. In Eq. Bll, 3^u, as 
defined in the main paper, is the Ising gauge flux through 



the plaquette of the direct lattice with the dual lattice 
point (a) at the centre. The implications of Eqs. BIO and 
Bll are straight forward. Eq. Bll states that flips the 
Ising flux and indeed is the vison creation operator. On 
the other hand, when combined with Eqs. B9, Eq. BIO 
shows that the VBS order parameter is indeed bilinear 
in the vison creation operators 

(Si . Sj) ^ {plpl). (B12) 

This completes our derivation for the VBS order param- 
eter. We also note that there are several alternative ways 
to derive the same result. An example is to start from 
the fermionic description of spinous as done by Senthil et 



Appendix C: The critical field theory 

The critical theory given by Eq. 33 where different 
terms are given by Eqs. 18,31 and 32 respectively. Now 
the flrst term in the vison action in Eq. 31 may be written 
as, 

TT 

Pab COS {4>a - 4>b) = COS (0a " 0& + -^(1 " Pab))- (CI) 

Within Villain approximation^^, the cosine term becomes 
(we henceforth use the notation A^^b = 4>a — i'b) 



oo 

COS (A0ai, + f ^ 



e 

mah — — 0O 

where mab is an integer defined on the links of the dual lattice. The partition function then written as 



^c^^E E E jT[d'^:5{\.,\'-l)\{dc^a e-'^-'^-'os (C3) 

{cr.j} {/lab} im^b} i a 

where, 

S„ = |^(A0a5 + |(l-Ma6))', (C4) 

{a,b) 

while §2 and §cs remain unchanged. Introducing the vison current field Jab on the links of the dual lattice through 
Hubbard-Stratonovich decoupling of (A0a6 + §(1 ^ Mab)) we get: 

^^=E E E [l[d'-iS{\z,\'-^)]ldc^a e-^-^-^- (C5) 

{<^,j} {f^ab} {Jab} i 

where, 

^v=Y.^+^J-b(^C^ab + \{l-Pab)) (C6) 

(ab) " 

Regrouping the terms in the critical action, Sp, we get: 

= + E # + 'J^^^ + f E(l - P<^^) 1^11°^) (C7) 

{ab) " {ab) ^ ^ 
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The sum over fiab can now be freely performed over each 
link of the dual lattice independently. This yields the 
constraint 



— Jab + 2naf, 



(where Uat is an integer) or, 



and §1, is given by: 



S.„ = ^ ^ + I Jab {A4>ab 



(ab) 



21 



(C8) 



(C9) 



(CIO) 



Finally integrating over the vison fields (/)o gives a further 
constraint on vison current fields Jab 



"^bJab = 



Appendix D: Breaking §U(2) to U(l) 



(Cll) 



To find the correct symmetry allowed term that breaks 
the second §U(2) to U(l) we consider the transformation 
of the X fields and the fields under various symmetries 
of the microscopic Hamiltonian. This may be obtained 
in several steps, due to the various transformations done 
to the microscopic Hamiltonian. A fruitful starting point 
is the transformations of the spinon fields z given by Eqs. 
9 and 10. We also use the same notation to denote the 
various transformations. 

The transformations of the spinous immediately im- 
ply the transformations for u. The transformation under 
spin rotation is given by 



(Dl) 



where F is a superposition of the generators of the first 
§U(2). For the rest of the transformations we have 



T 



' af3 

Va -i {Txt^y)af3 



Vl3 



(D2) 



Since the gauge fields are minimally coupled to the spinon 
fields their transformation and also the transformation of 



C fields are naturally implied by the transformation of the 
spinon fields. They arc given by 



"Ry 

7 



Up u,^ 

a-x ax, CLy -ay, 

ax — > —ax, ay — > —ay, a 
a^L ^ a-ti 



(D3) 



Combining all these, finally we have the transformation 
for the X fields given by: 

Ta : Xa^ (e-'^^A xp 

\ / ap 

^3,3 : Xa -> {TxfJ-y)apXp 

7 ■■ Xc.^{Tx)^pX} (D4) 

Thus to break the second §U(2) to U(l) (refer Eq. 52) 
we may consider a generic term (It is interesting to note 
that Ty generates translation for the x fields) 



(D5) 



where rt is a positive integer. Clearly such a term breaks 
the second §U(2) to U(l) and hence produces the correct 
symmetry of the microscopic Hamiltonian. To determine 
the lowest value of n consistent with various discrete sym- 
metries we consider the transformation of this term under 
various symmetries as before: 



Ta 

T 



ix* 
ix* 



ix* 
ix* 
x) ^ - ix* 



■ X) 
Ty ■ X) 



Ty ■ X) 
Ty ■ X) 



■X) 



(D6) 



Thus the linear term {n ~ 1) is forbidden by by time 
reversal symmetry. Thus the lowest order term allowed 
by all the microscopic symmetries is for n = 2, i.e. 



ix* ■ A2 • xf 



(D7) 
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